High efficiency switching using graphene based electron 'optics' 
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The absence of a band-gap in graphene hmits the gate modulation of its electron conductivity, 
both in regular graphene as well as in PN junctions, where electrostatic barriers prove transparent 
to Klein tunneling. We demonstrate a novel way to directly open a gate-tunable transmission gap 
across graphene PN junctions (GPNJ) by introducing an additional barrier in the middle that re- 
places Klein tunneling with regular tunneling, allowing us to electrostatically modulate the current 
by several orders of magnitude. The gap arises by angularly sorting electrons by their longitudinal 
energy and filtering out the hottest, normally incident electrons with the tunnel barrier, and the 
rest through total internal reflection. Using analytical and atomistic numerical studies of quantum 
transport, we show that the complete filtering of all incident electrons causes the GPNJ to act as a 
novel metamaterial with a unique gate-tunable transmission-gap that generates a sharp non-thermal 
switching of electrons. In fact, the transmission gap gradually diminishes to zero as we electrostat- 
ically reduce the voltage gradient across the junction towards the homogeneous doping limit. The 
resulting gate tunable metal-insulator transition enables the electrons to overcome the classic room 
temperature switching limit of ksTlnlO/q ^ 60 mV/decade for subthreshold conduction. 



Interest in graphene electronics is largely driven by its 
impressive material and electronic properties [1-3]. The 
two graphene bands derive from bonding-antibonding 
combinations of neighboring carbon Pz dimer ('pseudo- 
spin') basis sets belonging to the same two-dimensional 
crystallographic point group. The corresponding gapless 
low-energy excitations generate an ultralow electron ef- 
fective mass, while the orthogonality of the pseudo-spins 
suppresses 1-D back-scattering [^ , resulting in an incred- 
ibly high electron mobility [sl-Q • Unfortunately, the gap- 
lessness also makes its switching properties quite mod- 
est [8] , the subthreshold current changing linearly rather 
than exponentially with voltage. Conversely, opening a 
bandgap structurally by chemical modifications or quan- 
tum confinement [y-11] fundamentally reduces the mo- 
bility due to an asymptotic constraint on its short wave- 
length behavior il2|, rendering high efficiency switching 
in graphene a considerable challenge. 

The phase correlation between conduction and valence 
band states allows uniquely novel forms of electron flow 
in graphene. At a graphene p-n junction, the scattering 
of the individual pseudospins leads to electron trajecto- 
ries that are reminiscent of E-M wave scattering at di- 
electric interfaces. Since the radius of the 2-D Fermi sur- 
face varies with local potential across an interface, elec- 
trons with large angles of incidence at the higher electro- 
statically doped side are unable to conserve their trans- 
verse quasi- momentum components across the junction 
and thus undergo total internal reflection. Those within 
a critical angle 'refract' into the opposite side, diverging if 
the doping has the same sign (i.e. n+n or p+p junctions), 
or converging if opposite (pn junctions). A varying elec- 
trostatic doping thus turns graphitic electrons into quan- 
tum mechanical analogues of negative index metamate- 
rials [13|, the linear E-k making the electron trajectories 
non-dispersive over a finite energy window imposed by 



temperature and drain bias. Such electron 'optics' al- 
lows us to modulate the conductivity of a graphene sheet 
with split gates (l#J17l or build electronic analogues of 
fiber optic cables [18]. However, normally incident elec- 
trons well below the critical angle can directly tunnel 
between bands ('Klein tunneling') even for large voltage 
gradients, once again reducing the ON-OFF ratio. 

In this letter, we explore a graphene junction where 
electrons injected by a point source are spectrally sep- 
arated by a local gate, and those refracted across the 
junction are collected with an extended drain (Fig. 1). 
When electrons are injected from a higher electrostati- 
cally 'doped' to a lower doped side, the ones with high 
incident angle, low longitudinal energy are eliminated 
by total internal reflection. An additional tunnel bar- 
rier removes the normally incident, hotter electrons as 
well, resulting in a voltage- dependent transmission gap 
that allows significant gate modulation of the electron 
current (Figs. [2][3|). Furthermore, as we vary the gate 
voltage on the drain side, moving progressively towards 
a homogenous doping across the junction, we see an 
effective upconversion of the local voltage, leading to 
a subthreshold swing lower than the textbook thermal 
limit of kBT\nlO/q. The reduction in swing, that ulti- 
mately limits low energy switching, arises from the volt- 
age tunability of the transmission gap, and in fact col- 
lapses the bandgap in the homogeneous doping limit, re- 
quiring the 'valence band-edge' to slow down while the 
'conduction band-edge' to catch up. We thus have a 
unique non-thermal switching mechanism compared to 
other low subthreshold swing devices proposed J19l - [21| , 
arising specifically from a gate-tunable metal-insulator 
transition that bypasses the traditional limitations of 
tunnel based switches such as low ON current. 

Normal transmission in a GPNJ. The wavefunctions 
for graphitic electrons follow from the Dirac Hamilto- 
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FIG. 1. (a) A point source angularly spreads electrons while an extended drain collects those that refract around a tunnel 
barrier. Split gates control the two graphene segments separately. For the energy band diagram (b) with |Vgi|> |Vg2| and 
very small drain bias, (c) an atomistic current density plot calculated with the Non-Equilibrium Green's function (NEGF) 
formalism illustrates that current flows only within the critical angle, Oc = sin~'^{VG2/VGi) due to direct band-to-band (Klein) 
tunneling, (d) Schematic Snell's Law predictions as well as (e) fully atomistic quantum calculations illustrate how a barrier 
eliminates transmission within the critical angle when Oc < Ob, where Ob is the occlusion half angle created by the barrier. 
The non zero transmission occurs only for Oc > Ob where electrons can refract around the barrier (f). The elimination of 
transmission by coupling angular filtering with energy filtering ultimately generates the high switching efficiency of GPNJs 



nian, H{k) = hvp^-k^ where a = {(Jx^cFy) are the PauH 
matrices and vp ^ l^^ms~^ is the Fermi velocity for elec- 
trons in graphene. The wave- functions consist of a plane 
wave part with Fermi wavevector magnitudes kFi,F2 ori- 
ented along incident (^i), reflected (^i) and transmitted 
angles (^2) on the two sides, while their Bloch-like atomic 
parts act as two-component spinors whose phases are 
determined by the angles and overall signs determined 
by the band-indices. Matching the plane wave phases 
across the boundary, we get the corresponding Snell's 
Law kpisinOi = —kF2sin02 with a refractive index ratio 
n = —kF2/kFi' Matching the spinor (Bloch) components 
thereafter gives us the transmission probability 



T{ei,E,VG) = ei0c-ei)- 



cos6icos92 
91+92 



(1) 



where 9 is the unit step function, O2 is related to Oi 
by Snell's law, Vg={Vgi-, Vg2} is a particular set of gate 
voltages and Oc = sin~^{n) is the critical angle. The 
equation gives zero transmission (total internal reflec- 
tion) for Oi > Oc^ unit transmission for a homogenous 
sheet {Oi = ^2), and focusing with T{Oi) = cos^Oi for 
symmetric GPNJs (^1 = —O2)' The transmission can also 
be plotted vs. longitud inal energy El for a given total 
energy E = ^J E\ + E\. The critical angle Oc eliminates 
the obliquely incident electrons with high transverse and 
low longitudinal energy, so that T{El) looks like a high- 
pass filter^ transmitting only states with El > E cos Oc- 

The energy dependence of T enters through the terms 
related by Snell's law, while the gate voltage depen- 
dence enters through the positions of the charge neutral- 



ity points on both sides of the junction. The transmission 
here is that of a single mode {kx^ky} with incident angle 
Oi = tan~^{ky/kx) and energy E = hvF\/{k'^ + ^^) = 
\/ E'j^ + E'^. The average transmission of all modes at a 
given total energy E is found by integrating T from Eq. 
[TJover all incident angles and then normalizing (Fig. [2^). 
As E approaches the neutrality point on the refracted 
side, more and more modes suffer total internal reflection, 
until at —qVG2 (0.4 eV in this case), the critical angle van- 
ishes and all modes are reflected. Since this happens at 
a single voltage where all incident modes align with the 
Dirac point on the refracted side, the ON-OFF ratio is 
quite poor and in fact comparable to a regular graphene 
transistor. In effect, electrons are injected over a window 
of energy set by the drain bias, creating a finite cone over 
which transmission will take place (Fig. [It) and degrad- 
ing the ON-OFF. Since electrons with higher longitudinal 
energy see a lower barrier, eliminating transmission from 
the high longitudinal energy electrons is critical to higher 
efficiency of the Klein tunnel switch. 

Creating a transmission gap. Fig. [1^ shows the 
schematic diagram of the device we propose. The source 
injects electrons at all angles. If we pattern a barrier 
at the center of the graphene sheet, the normally inci- 
dent electrons will be eliminated, so that the transmis- 
sion probability with barrier can be written as 



Tb{9i,E,Vg) = e{9i-eB)T{euE,VG) 



(2) 



where Ob = tan~^D/2L is the occlusion half angle sub- 
tended at the source by the barrier, D is length of the 
rectangular barrier and L is its perpendicular distance 
from the source. Atomistic NEGF calculations in Fig. [1^ 
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FIG. 2. (a) GPNJ Transmission vs. Fermi energy from Eq.[T] 
for Vgi=+0.4V and Vg2=-0.4V. At all energies, except when 
aligned with the Dirac point on the p side, the zero angle mode 
transmits with unit probability, making the average transmis- 
sion nonzero, (b) The circles represent constant Fermi energy 
slices of the Dirac cones across the junction, corresponding to 
points on the T(E) curve. Arrows indicate the maximum ky 
values that transmit across the junction, and thus define the 
critical angle, (c) Modification of the transmission spectrum 
from part (a) after putting the barrier; transmission within 
the yellow region is forbidden due to barrier, (d) Replotting 
T vs longitudinal energy El at specific total energies E, indi- 
cating an effective band-pass behavior for states outside the 
T(E) transmission gap. 



confirm how modes near normal incidence are reflected 
by the barrier, while those incident at a large angle are 
eliminated by total internal reflection. The transmission 
now acts as a band-pass filter along the longitudinal en- 
ergy axis El (Fig. [2]i), allowing only electrons with an- 
gles Ob < Oi < Oc^ in other words, within a longitudinal 
energy window E cos Oc < El < E cos Ob to transmit. 
Since no electrons can transmit when Oc < Ob, there is 
a range of energies for which we get a transmission gap. 
The gate voltage dependance of Oc leads thereby to a 
voltage-tunable transmission gap. 

Gating the transmission gap: a continuous metal- 
insulator transition. We have shown that filtering of elec- 
trons with high longitudinal energy leads to a transmis- 
sion gap. One can extract an effective 'valence band- 
edge' Ey when this happens, given by the condition 
sinOB = sinOc = —{Ey + qVG2)/{Ev + qVgi)- One can 
get a similar expression for the 'conduction band-edge' 
Ec, and thus the effective transmission gap Eq 



Ev,c = -qVG2 T q 



(Vbi - Vb2) sin i9b 
1 ± siuOb 



EG = Ec-Ev=q{VGi-VG2)^^^ (3) 

CO 8"^ Ob 

Note that there are two distinct contributions to the re- 
sulting transmission plot (Fig. [3^)- (i) an overall shift 



given by —qVG2, and (ii) a voltage modulation of the 
band-gap Eg- The first term will lead to the electrostatic 
gating effect seen for regular band-gapped semiconduc- 
tors (for good gate control geometries, this amounts to 
kBT In 10 /q ^ 60 mV/decade), while the latter term will 
give a deviation from this textbook result, effectively cap- 
tured through a gate tunable, continuous metal-insulator 
transition associated with voltage dependence of the 
transmission gap Eg- As expected, the band-gap van- 
ishes in the absence of either a tunnel barrier {Ob = 0), 
or a voltage gradient across the junction {Vg2 = Vgi, 
^c = 900). 
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FIG. 3. (a) Transmission vs total energy for different gate 
voltages. The transmission shifts as we vary the voltage Vg2, 
but is accompanied with a change in transmission gap; (b) 
I-Vg2 of the device, showing a sharp increase in current mod- 
ulation with barrier, and a subthreshold swing that is larger 
than the textbook limit at the PN end and smaller at the 
N+N end. Here, Vgi =+lvolt, Ob = 20° and a drain voltage 
of 0.2V. 

Fig. [3] shows the current extracted by integrating the 
analytical transmission formula (Eq. [2]) using the Lan- 
dauer formahsm for various gate voltages Vg2 while Vgi 
is kept fixed. While the ON current reduces a little, the 
tunnel barrier primarily reduces the OFF current (Fig. 
4), dropping it by nearly five orders of magnitude. Over- 
all, an ON-OFF ratio of ^ 10^ is achievable with just 
0.25 V change in gate voltage on the NN regime with Ob 
of20^ 

The unusual transconductance and subthreshold swing 
in Fig. [3] can attributed to the fact that transmission 
gap changes everytime the gate voltage Vg2 is changed. 
Eq|3] gives us the edges of the transmission gap, so 
dEcy /dVG2 = (^gQ/[^ T siuOb], where we have now in- 
troduced the capacitive gate transfer factor ag. In other 
words, the subthreshold swing, ultimately dictating the 
efficiency of low voltage switching, will be 

S = {dlogioI/dVG2)-' = M:1^[i ^ sinOB] (4) 

agQ 

which is less than 60mV/dec ('-' sign above) at the NN 
end, and more than 60mV/dec ('+' sign) at the NP end. 
At heart of the unusual behavior is the coupling of angu- 
lar and energy filtering in the barrier driven Klein tunnel- 
switch, which differs from most materials in that its ef- 



fective band-gap can be collapsed with a gate voltage. 

The tunability of transmission gap is reproduced with 
a fuhy atomistic NEGF tight-binding calculation of cur- 
rent flow across a GPNJ sheet involving ~ 100,000 atoms 
(Figs. 1, 4). The transmission gap is slightly lower than 
the analytical predictions owing to the presence of leak- 
age currents created by edge reflection of the electron 
waves. Such effects are expected to be substantially re- 
duced if we edge passivate large sheets of graphene. In- 
deed, selectively imposing absorbing self-energy matrices 
minimizes such standing waves, causing the computed 
transmission gap to approach our analytical predictions. 

The physics of regular tunneling coupled with Klein 
tunneling opens up the opportunity for high performance 
low power switching based on graphene. Our results are 
predicated upon Snell's Law, in other words, the preser- 
vation of the transverse quasi- momentum across the junc- 
tion. Various non-idealities can compromise these ef- 
fects, such as width- independent leakage currents arising 
from edge reflections and diffraction effects at the barrier 
edges. These non- idealities are already accounted for by 
our numerical calculations and dominate primarily for 



long wavelength electrons that are, in fact, eliminated 
by the tunnel barrier. For a l/iTn x l/iTn sheet with an 
occlusion half angle Ob = 20^, a ~ 180 nm x Inm rectan- 
gular void generates a tunneling probability that is con- 
siderably less than 10~^. Experimental success depends 
on the mobilities of graphene samples, and the quality 
of the tunnel barrier created either using electron beam 
lithography or by directed assembly of insulating molec- 
ular species on graphene, as demonstrated by Hersam 
et al [23|. Further numerical simulations are needed to 
study the robustness of these effects in the face of charge 
puddles that could misalign the Dirac points within each 
graphene segment, finite size effects at the point con- 
tact, line-edge roughness at the electrodes, finite junction 
widths imposed by the Debye lengths, and pseudospin 
mixing by edge states with various degrees of chemical 
passivation. 
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